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摘 要 本 研究 基于 二 氧化 碳 气 体 v3 基础 吸收 带 的 直接 吸收 光谱 测量 乙烯 同 轴 扩 散 火 焰 的 温度 场 。 实 验 采用 波长 
为 4.17 微米 的 室温 带 间 级 联 激光 器 (interband cascade laser), 测量 此 吸收 带 多 个 吸收 峰 通 过 轴 对 称 火焰 的 路 径 积 分 吸 
收 ， 并 通过 阿 贝 尔 道 变换 方法 ， 反 演 得 出 火焰 不 同 高 度 吸 收 率 沿 火焰 径 向 的 分 布 ， 进 而 参考 HITRAN 数 据 库 通过 最 
小 二 乘法 拟 合 实验 数据 得 到 火焰 中 二 氧化 碳 的 温度 场 分 布 。 为 了 提高 阿 贝尔 转化 的 精度 ， 本 文采 用 正则 化 阿 贝尔 转 
化 方法 ， 并 通过 数值 模拟 验证 了 正则 化 方法 的 可 行 性 。 
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Temperature Distribution of Axisymmetric Diffusion Flame Measured by 


Mid-infrared laser 


HUANG Yan ZHANGGuo-Yong LIU Xun-Chen QI-Fei 
(School of Mechanical Engineering, Shanghai Jiao Tong University, Shanghai, 200240, China) 
Abstract: Temperature distribution of axisymmetric ethylene-air coflow flame was measured by direct 


absorption spectroscopy at the v3 band of carbon dioxide. A room temperature interband cascade laser 
(ICL) was used to continuously scan wavelength at the bandhead region near 4.17 um. Along the 
axisymmetric coflow diffusion flame, the path integrated absorption signal was inversed to linear 
absorption coefficients through Abel inversion, which was then non-linear least square fitted with 
HITRAN database to reconstruct the temperature and carbon dioxide distributions in the flame. In order to 
increase the accuracy of Abel inversion, Tikhonov regularization method was used and its feasibility was 
verified through numerical simulation. 
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0 引言 


层 流 同 轴 扩 散 火 焰 作 为 一 种 标准 火焰 ， 常 用 于 
碳 烟 生 成 和 燃烧 化 学 的 研究 5 。 火 焰 的 温度 场 和 
燃烧 产物 浓度 信息 是 重要 的 火焰 参数 ， 决 定 了 碳 烟 
的 生成 、 氧 化 和 燃烧 反应 的 速率 。 热 电 侦 常 被 用 于 
气体 温度 的 测量 。 然 而 ,对 火焰 采用 热电 偶 测 温 时 ， 
热电 偶 焊 点 周围 的 辐射 热传导 可 导致 焊 点 与 周围 气 
体 之 间 存 在 热 不 平衡 ， 且 用 于 制作 高 温 热 电 偶 的 材 
料 在 火焰 环境 中 的 不 稳定 会 增加 辐射 传 热 校正 中 相 
关 的 不 确定 性 中。 可 调谐 二 极 管 激 光 吸 收 光 谱 
(TDLAS) 技 术 可 实现 火焰 中 温度 和 气体 浓度 的 原 位 
测量 9。 火焰 中 的 二 氧化 碳 在 4.17 微米 的 中 红外 
波段 有 较 强 的 吸收 。 在 该 波长 附近 ， 无 需 考 虑 空气 
中 其 他 气体 物质 或 火焰 中 的 其 他 气体 燃烧 产物 对 激 
光 的 吸收 干扰 ， 因 而 便于 测量 且 精 度 较 高 54。 直 接 
吸收 光谱 技术 测量 的 是 沿 激光 路 径 的 积分 信号 ， 对 
于 非 均 匀 分 布 的 火焰 参数 ， 通 常 需要 采用 层 析 方法 
对 激光 路 径 上 各 个 位 置 处 的 火焰 参数 进行 重 构 门 。 
对 于 轴 对 称 火焰 ， 如 层 流 平面 预 混 火焰 入 和 层 流 同 
轴 扩 散 火 焰 ， 温 度 和 气体 浓度 等 火焰 参数 呈 轴 对 称 
径 向 分 布 ， 因 此 可 通过 测量 距离 火焰 中 心 不 同 距离 
的 路 径 积 分 信号 和 阿 贝 尔 逆 变 换 方法 有 反 演 重 构 。 

阿 贝 尔 道 变换 在 科学 和 工程 上 有 具有 广泛 应 用 ， 
目前 已 发 展 了 许多 不 同 的 逆 变 换 方 法 ， 其 中 最 常用 
的 是 Onion-Peeling( 下 文中 称 OP) 和 Abel Three 
Point( 下 文中 称 ATP) 方 法 名。 然而 此 类 方法 由 于 其 
相应 矩阵 方程 的 病态 特性 由， 实验 信号 的 微小 误 
差 在 反 演 求解 过 程 中 会 被 放大 ， 使 得 反 演 结果 明显 
偏离 真实 值 。Tikhonov 正则 化 方法 被 证 实 可 以 有 效 
地 减少 阿 贝尔 道 变换 的 误差 &0。 正 确 使 用 
Tikhonov 正则 化 方法 的 关键 在 于 确定 合适 的 正则 化 
参数 ， 其 中 L-curve 方法 是 常用 的 确定 正则 化 参数 
Kae, 

使 用 TDLAS 技术 结合 阿 贝 尔 逆 变换 方法 重建 
轴 对 称 火 焰 温 度 场 的 工作 多 针对 于 平面 火焰 8, 
还 未 有 使 用 该 方法 对 扩散 火焰 温度 场 的 研究 。 本 文 
通过 数值 模拟 验证 了 正则 化 方法 的 适用 性 ， 并 通过 
4.17 微米 波长 的 室温 带 间 级 联 激光 (ICL) 测 量 乙烯 
同 轴 扩 散 火 焰 中 二 氧化 碳 的 路 径 积 分 吸收 信号 ， 结 
合 阿 贝尔 道 变换 和 HITRAN 数据 库 拟 合 获得 了 乙 
烯 扩散 火焰 的 温度 分 布 。 


1 基本 原理 和 数值 模拟 
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1.1 直接 吸收 光谱 原理 

TDLAS 技术 建立 在 比尔 郎 伯 吸收 定律 的 基础 之 
上 ， 当 频率 为 的 激光 束 穿 过 长 度 为 L 的 被 测 气体 
时 ， 激 光 光 强 的 透 过 率 开 可 由 沿 光 程 的 吸收 率 积 4 
得 到 : 


T, = (2), =expC- 人 cl 加 () 


all, v= XD):S TO) fv ),TO) (2) 


式 中 ，7 和 7 分 别 为 入 射 激光 强度 和 透 过 激光 
WE, (D Al 位 置 处 的 被 测 气体 的 摩尔 分 数 ， 
a(l, v) 为 ! 位 置 处 的 被 测 气体 对 频率 为 的 激光 的 
吸收 率 ， 吸 收 系数 S, 是 跃迁 的 谱 线 强度 ，/, 是 跃 
E j 的 线 型 函数 ，v ;是 跃迁 j 的 中 心 跃迁 频率 。 线 
强 S, 与 温度 7 的 关系 可 由 下 式 确定 : 

E ODEN 
I Oar T ETT 3) 
1—exp(—Acv,; / kT) 
1-exp(-hcv, , / kT) 

AF, h 为 普 朗 克 常 数 ，c 为 光速 ,，k WEKA 
BWM, N 为 参考 温度 (通常 为 296 K) E, 为 吸收 
跃迁 的 低能 级 能 量 ，Q(D 是 配 分 函数 “了 ]。 从 上 式 可 
以 看 出 ， 线 强 仅 与 温度 有 关 。 

本 研究 中 的 常 压 扩散 火焰 测量 的 光谱 线 型 可 通 
过 Voigt 线 型 来 描述 。Voigt 线 型 由 高 斯 线 型 和 洛 伦 
效 线 型 的 卷 积 得 到 ， 可 由 下 面 的 公式 得 到 : 


Vu,o,7)=| Gu ,o)L(u-u',y)du (4) 
其 中 , u 为 相对 频率 (实际 频率 -中 心 频率 )，a 和 yy 分 
别 是 高 斯 线 型 和 洛 伦 效 线 型 的 半 高 半 宽 (HEWHM) , 高 
斯 线 型 和 洛 伦 兹 线 型 可 分 别 由 下 式 给 出 : 
Guon — w 2 (5) 
TO oO 
7 

L(y,7)=—— 6 

(u,y) ETS (6) 


AWRA, EARE EREE Oo 和 yy 可 由 下 式 
给 出 : 


ee Vy [2 In2N ,kT (7) 
C M 


T Noir 
y=P- z(=) (8) 


KP, N, 为 阿 伏 加 德 罗 常 数 ，M 为 分 子 摩尔 质量 ， 
/是 常温 常 压 下 的 碰撞 展 宽 的 半 高 半 宽 ，n, 与 温度 
相关 的 指数 系数 。HITRAN 数据 库 提供 了 不 同 跃迁 
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WEF ya Mn, AVRO. KEREKE 
下 进行 ， 因 此 吸收 率 主要 取决 于 火焰 温度 和 被 测 二 
氧化 碳 气 体 的 浓度 。 


12 阿 贝 尔 逆 变换 方法 
图 1 是 轴 对 称 火焰 某 个 高 度 的 水 平 剖 视图 ， 根 
据 比 尔 即 伯 定律 和 阿 贝尔 积分 方程 ， 对 于 某 一 频率 
vV， 可 以 得 到 吸收 率 alr, v) 的 路 径 积 分 测量 值 
P(x,v) 与 其 分 布 之 间 的 关系 : 
A e en (9) 
vV<v<v,, —R<x<R 
AF y 和 vw 分 别 是 连续 激光 扫 面 频率 的 下 限 和 上 
限 。 阿 贝尔 方程 10) 具 有 解析 解 : 
1 rR[dP(x,v)/ dx 
a(ry)=-—[ aE (10) 
vy<v<v,, —R<r<R 
RP dP(x,v)/ dx 是 测量 路 径 积分 值 P(x,v) E x h 
方向 上 的 导数 。 然 而 在 实际 测量 中 ， 通 过 有 限 差分 
方法 获得 P(x,v) 的 导数 来 计算 alr, v) 常常 会 放大 
通过 式 (10) 直 接 求解 的 方法 不 具有 实 


x 


图 1 轴 对 称 扩散 火焰 激光 测量 


Fig.1 Laser measurement of axisymmetric coflow Flame 


许多 替代 的 阿 贝尔 算法 相继 被 提出 。OP 方法 将 
火焰 水 平 剖面 划分 为 N 个 宽度 相等 的 圆 环 ， 并 假设 
在 每 个 圆 环 内 火焰 参数 相等 ， 将 阿 贝 尔 逆 变 换 问题 
转化 为 求解 线性 方程 组 
Apa =P (11) 
HH: @={),%,0,,...,Ay,} 表示 从 圆心 向 外 ， 各 
个 圆 环 内 的 吸收 率 的 值 ，P= (P,P, P,,..., Py} 表示 
距离 圆心 不 同 的 距离 的 路 径 积 分 信号 ，4o 是 尺寸 为 
NxN 的 几何 关系 矩阵 。ATP 方法 转化 精度 更 高 ， 其 
考虑 相 邻 路 径 的 偏 导 门 , 并 将 阿 贝尔 逆 变 换 转 化 为 求 
解 类 似 的 线性 方程 组 外 
Air æ =P (12) 
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1.2.1 ”病态 和 矩阵 问题 

由 于 OP 和 ATP 方法 所 构造 矩阵 的 病态 特性 
实验 测量 的 路 径 积 分 信号 的 微小 噪声 会 导致 求解 值 
与 真实 值 偏 差 较 大 。 在 求解 方程 Aa = PIN, ARE 
阵 4 进行 奇异 值 分 解 趾 ， 可 以 得 到 

4=UZ =) uoy; (13) 


` 


AP: UR V HAG N 个 正 交 列 向 量 u, Flv,» È 
是 包含 矩阵 4 WaT HELO, MIME. ATHERE U 
和 严 的 正 交 性 ， 方 程 的 解 可 以 写成 

= hu p x u; Po , , MOP 14 
a pe: pa aes v,+ 7 v,) (14) 


AP: 万 是 由 实验 测量 获得 的 路 径 积分 吸收 信号， 
Dy 是 没有 实验 误差 时 理想 的 路 径 积分 吸收 信号， 
Op = p 一 po 是 实验 测量 结果 的 误差 。 从 (14) 式 可 以 
AH, RE A 的 小 奇异 值 会 放大 实验 误差 的 影响 ， 
从 而 降低 阿 贝 尔 转 化 的 精确 度 。 
1.2.2 Tikhonov 正则 化 方法 

对 于 稳定 的 火焰 ， 其 温度 场 和 气体 浓度 分 布 应 
当 是 连续 平滑 的 。 一 种 有 效 的 正则 化 方法 是 对 上 述 
和 矩阵 方程 添加 约束 项 ， 即 添加 平滑 约束 条 件 ， 使 用 
修正 的 方程 进行 阿 贝 尔 逆 变 换 求解 "由 。 因 此 ， 可 
将 原来 的 矩阵 方程 4c = 已 修正 成 以 下 形式 : 


E o 


式 中 : 工 是 用 于 约束 平滑 条 件 的 拉 普 拉 斯 矩阵 ,4 是 
控制 正则 化 程度 的 正则 化 系数 。 和 矩阵 工 的 形式 如 下 : 


1 -1 0 … 0 

Le 0 1 -l - (16) 
: = 
0 … 0 1 -l 


EREL HTE 中 的 相 邻 项 平滑 过 渡 , EE L 
的 行 数 k 可 根据 不 同 阿 贝尔 转化 的 对 象 以 及 阿 贝 尔 
转化 矩阵 的 不 同 而 做 调整 ，k=1,2...N 一 1 XF 
OP 方法 , k=1 (对 方程 加 约束 项 A Q0-01)=0 时 即 
可 以 达到 较 好 的 转化 效果 ; 而 对 于 ATP 方法 ,Kk 的 
选择 对 正则 化 效果 的 影响 较 小 。 下 文中 模拟 中 ， 对 
OP 方法 ， 取 k=1 ， 对 ATP 方法 ， 取 k=N 一 1 。 
由 于 方程 (15) 是 超 定 方程 , 可 通过 最 小 二 乘法 求 
最 优 解 ， 即 求解 下 列 方程 : 


| 4" ary i, ele pa (17) 
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简化 成 
| 44+ VE L]a=A'P (18) 


NA 


max , (Kk, ) 


> 


Pa 
图 2 用 于 确定 正则 化 系数 的 L-Curve 

Fig.2 L-Curve method for selecting regularization parameter 

采用 正则 化 的 阿 贝 尔 转换 方法 时 ， 正 则 化 系数 
4 的 选择 十 分 关键 : 4 过 小 ， 正 则 化 效果 不 明显 ; 
而 4 过 大 时 ， 尽 管 结果 曲线 光滑 ， 但 准确 性 无 法 保 
证 。L-curve 方法 被 证 实 是 比较 简单 有 效 的 确定 正则 
化 参数 的 方法 中。 图 2 表示 了 L-curve 方法 的 基本 
原理 , 对 于 不 同 的 4 的 值 , 绘 出 对 应 不 同 求解 值 x; 的 
BAL ty Ay log ao (IIAx, 一 8lls) 纵 坐标 为 logso(llExa ll) 
的 各 点 连 成 的 曲线 。 FETE [lars lllr 一 如 :分 别 是 
Lx Ar,- DHLA WTE A, mR lat 
大 ， 则 该 逆 变 换 仍 然 对 测量 噪声 较为 敏感 ， 如 果 
Ax, 一 8 上 :过 大 ， 则 逆 变 换 的 结果 精度 较 低 。 如 图 2 
所 示 , L-curve 出 现 最 大 曲率 时 , lix || mlr — bll 
都 较 小 ， 是 较为 理想 的 正则 化 系数 。 

在 本 实验 中 ， 需 要 对 coflow 火焰 的 多 个 高 度 处 
的 测量 数据 进行 阿 贝 尔 转化 ， 正 则 化 系数 的 自动 化 
选择 可 以 简化 数据 处 理 过 程 。 通 过 有 限 差分 求 导 可 
以 计算 曲线 曲率 ， 找 到 曲率 的 最 大 值 对 应 的 正则 化 
系数 ， 从 而 实现 正则 化 阿 贝尔 逆 变 换 的 自动 化 运算 : 
k, = pall bath 7 as 

OROI 

式 中 :pa = logs (llAx_ — bll), na = log oti ll). 
1.2.3 ”数值 模拟 验证 


在 coflow 火焰 的 较 低位 置 处 ， 径 向 吸收 率 主要 
呈 双 峰 分 布 ， 而 在 火焰 的 较 高 位 置 处 ， 径 向 吸收 率 
主要 呈 单 峰 分 布 。 为 了 便于 进行 模拟 计算 ， 我 们 构 
造 了 如 下 图 3 所 示 的 两 种 分 布 : 


1 2/2 
A) =~ Us | (20) 
flr) = (1 sr) (21) 
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式 中 ，fi(7) 和 (7) 分别 为 假设 的 沿 火 焰 水 平 径 向 
分 布 的 单 峰 和 双 峰 分 布 函数 , 其 中 分 布 半径 x 作 了 归 
一 化 处 理 。 根 据 式 (9) 使 用 分 部 积分 方法 可 以 计算 出 
阿 贝 尔 路 径 积分 函数 : 


RŽ)" (22) 


4 2 \3/2 8 2\5/2 

Pœ =z x ) sl! x ) (23) 

RP, P(x) M P (x) 分 别 是 单 峰 和 双 峰 分 布 在 
距离 火焰 中 心 x 位置 处 的 阿 贝 尔 路 径 积 分 值 。 本 实验 
中 ， 考 虑 对 20 条 等 间距 的 路 径 积 分 值 进行 逆 变 换 求 
解 。 当 不 考虑 路 径 积分 值 的 实验 误差 时 ， 即 使 是 最 
简单 的 OP 方法 也 能 获得 较为 精确 的 转化 结果 ,如 下 
图 3 所 示 。 

模拟 实验 对 积分 值 加 幅 值 为 2% 随 机 噪声 ， 再 对 
加 噪声 的 积分 值 分 别 使 用 OP、ATP、 正 则 化 的 OP( 下 
文中 称 Tik-OP) 和 正则 化 的 ATP 方法 (下 文中 称 
Tik-ATP) 进 行 逆 变换 求解 并 与 原始 分 布 式 比较 , 验证 
不 同 逆 变换 方法 的 求解 精度 ， 典 型 的 转化 结果 如 下 
图 4 所 示 。 


fin) 


— exact 
0.02 e 。 Onion Peeling 


0.080 0.2 0.4 0.6 0.8 T.0 


r 
图 3 不 考虑 噪声 时 候 的 Abel 转化 


Fig.3 Abel inversion without noise 
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图 4 考虑 噪声 时 的 阿 贝尔 转化 结果 


Fig.4 Abel inversion with noise 


为 了 定量 分 析 误 差 ， 重复 进行 了 100 次 含 噪声 
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数值 实验 ， 利 用 100 次 平均 的 均 方 根 误差 来 比较 不 
同方 法 在 有 噪声 情况 下 的 转化 精度 : 


至 | Siw — Jia ll, 24 
MN #max( Soa) a 
—_ | 
=r 25 
E ins 100 + (Ems J ( ) 


HH, Em 为 均 方 根 误差 ，N 为 fion Sraa 中 的 项 数 ， 
这 里 取 N=20 ， 广 为 径 向 分 布 的 原始 设 定 值 ，/ 为 
阿 贝 尔 逆 变换 的 求解 值 , BY 100 次 模拟 求解 的 均 
方 根 误差 的 平均 值 ， 结 果 如 下 表 1 所 示 。 

表 1 均 方 根 误 差 平均 值 比较 


Table 1 Comparison of Mean Square Error 


Abel Methods Single Peak Double Peak 
OP 0.08738 0.04345 
ATP 0.03030 0.02187 
Tik-OP 0.00491 0.01527 
Tik-ATP 0.00520 0.01877 


从 表 1 可 以 看 出 ， 各 种 方法 按照 求解 精度 排序 
依次 为 : Tik-OP>Tik-ATP>ATP>OP 。 对 于 尺寸 为 
20x20 的 OP HERE, 本 文采 用 的 正则 化 方法 可 以 有 效 
提高 转换 精度 ， 而 对 于 ATP 矩阵， 正则 化 方法 对 精 
度 的 提高 有 限 。 

基于 上 述 模拟 结果 ， 在 处 理 整 个 扩散 火焰 的 实 
丛 数 据 时 候 ,可 采用 Tik-OP 方法 进行 阿 贝尔 逆 变 换 。 


2 实验 测量 


图 5 是 实验 装置 的 示意 图 。 本 实验 选用 波长 在 
4.17 微米 附近 的 ICL 激光 器 ， 通过 函数 发 生 器 产生 
的 正 驼 波 扫描 波 数 范围 在 2396.9 cm’! 至 2397.4 cm’! 
之 间 的 中 红外 激光 波长 。 激 光 器 产生 的 激光 通过 两 
个 分 束 镜 分 成 三 路 光 ， 即 分 别 由 三 个 探测 器 测量 信 
号 的 标准 具 波 长 校准 光路 、 参考 光路 和 coflow 火焰 
测量 光路 。 为 了 防止 激光 经 过 高 温 火 焰 区 域 发生 偏 
折 而 无 法 全 部 被 探测 器 接收 ， 在 探测 器 前 设置 聚焦 


镜 。 


压缩 空 
分 束 镜 C)H,/N> 
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图 5 实验 装置 示意 图 
Fig.5 Schematic view of the experimental apparatus 
实验 中 ， 步 进 电 机 控制 燃烧 炉 在 水 平方 向 和 坚 
直方 向 上 移动 ， 其 中 水 平方 向 测量 火焰 中 心 轴线 到 
火焰 外 缘 直 至 没有 吸收 信号 的 位 置 ,移动 间距 为 0.05 
em, 竖 直 方向 从 火焰 底部 到 可 见 火焰 尖 部 以 上 位 置 ， 
移动 间距 0.1 cm。 探 测 器 2 测 得 的 参考 光路 的 信号 
用 作 背 景 信号 1, ， 由 于 参考 光路 和 经 过 火焰 的 光路 
初始 激光 强度 不 同 ， 基 线 思 通过 除 以 距离 火焰 较 远 
位 置 处 的 两 路 信号 的 峰 峰 值 的 比值 进行 校正 。 经 过 
标准 具 的 激光 信号 具有 相 邻 的 峰 间隔 相同 波 数 的 性 
质 ， 结 合 HITRAN 数据 库 中 已 知 的 基本 固定 的 吸收 
峰 波 数 ， 可 标定 在 每 个 扫描 周期 内 的 激光 波长 信息 。 
图 6 示 出 了 火焰 中 某 个 测量 位 置 的 半 个 正弦 扫描 周 
期 内 的 原始 数据 ， 其 中 背景 信号 经 过 如 上 上 所 述 的 校 
正 ,。 图 6 中 可 以 看 出 在 火焰 温度 下 , 二氧化碳 对 4.17 
微米 的 中 红外 激光 的 吸收 较 强 。 
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图 6 实验 测量 的 信号 
Fig.6 Signals obtained in experiment 


为 比较 上 述 重 构 方法 对 单 峰 和 双 峰 分 布 的 实验 
测量 值 进 行 转换 的 准确 性 ， 分 别 对 实验 得 到 的 火焰 
高 度 Y=50 mm 和 Y=20 mm 位 置 处 的 吸收 信和 号 路 径 
积分 值 P(x) =—In(/ 1), 使 用 前 述 阿 贝尔 转化 方 
法 进行 逆 变 换 求解 并 进行 比较 ， 结 果 如 下 图 7 所 示 。 
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图 7 实验 值 的 阿 贝尔 转化 


Fig.7 Abel inversion of experimental data 


参考 模拟 验证 部 分 的 结论 并 结合 图 7 所 示 的 结 
果 , 可 以 判断 正则 化 的 OP 方法 可 以 获得 较为 光滑 准 
确 的 结果 。 因 此 对 火焰 的 每 一 测量 高 度 、 各 个 激光 
波 数 得 到 的 水 平方 向 路 径 积分 值 采用 Tik-OP 方法 ; 
行 阿 贝 尔 逆 变换 求解 ， 从 而 得 到 火焰 各 个 高 度 的 径 
向 吸收 率 分 布 ， 并 使 用 最 小 二 乘法 将 转化 后 的 吸收 
率 与 HTRAN 数据 库 中 的 数据 进行 拟 合 ， 从 而 得 到 
火焰 的 温度 分 布 。 

实验 中 选取 与 文献 2] 中 实验 相似 的 工 况 ， 混 合 
燃料 气体 中 乙烯 占 60%、 和 氮气 占 40%， 燃 料 出 口 流 
速 为 5 cms。 图 8-10 对 本 实验 的 结果 与 文献 中 中 公 
布 的 结果 进行 了 比较 ,其 中 火焰 外 围 较 远 处 的 部 分 点 
被 忽略 。 
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图 8 火焰 竖 直 剖面 的 二 维 温度 分 布 


Fig.8 Temperature distribution of vertical central section of flame 
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图 9 火焰 不 同 高 度 处 的 径 向 温度 分 布 
Fig.9 Radial distribution of flame temperature at several heights 
从 图 8-10 可 以 看 出 , 本 实验 测 得 的 火焰 温度 与 
文献 请 公布 的 结果 较为 相似 ;火焰 高 度 >30 mm 时 ， 
本 实验 获得 的 乙烯 火焰 温度 分 布 与 文献 站 中 结果 基 
本 一 致 ， 且 本 实验 的 温度 分 布 较为 平滑 。Y=25 mm 
时 , 本 实验 的 温度 与 文献 中 中 结果 相差 100 K 左右 ， 
但 是 径 向 温度 变化 趋势 较为 相似 。 本 实验 和 文献 站 
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中 的 实验 并 不 能 保证 完全 相同 的 实验 工 况 ， 例 如 燃 
烧 炉 的 尺寸 、 压 缩 空 气 的 流速 等 存在 差异 ， 可 能 是 
产生 不 同 测量 结果 的 原因 。 
3 结论 

本 研究 使 用 TDLAS 技术 并 结合 阿 贝 尔 转化 方 
法 重建 了 乙烯 轴 对 称 扩散 火焰 的 温度 场 。 为 提高 数 
据 处 理 的 精度 ， 采 用 了 Tikhonov 正则 化 方法 对 阿 贝 
尔 转化 矩阵 (OP 和 ATP 甜 阵 ) 进 行 约束 ， 并 进行 了 数 
值 模拟 比较 。 此 外 ， 为 了 验证 实验 测量 与 数据 分 析 
的 可 靠 性 ， 对 与 文献 [2] 中 实验 相似 的 工 况 进行 了 实 
验 ， 并 与 文献 中 结果 进行 了 比较 。 
DTikhonov 正则 化 方法 应 用 于 本 实验 中 尺寸 为 20 
的 OP 矩阵 时 ， 能 够 有 效 减少 路 径 积 分 噪声 对 转化 
结果 的 影响 ， 显 著 提 高 阿 贝尔 逆 变 换 的 精度 ; 
2) 测 得 的 乙烯 同 轴 扩散 火焰 的 最 高 温度 为 约 1850 
KK。 本 实验 获得 火焰 温度 分 布 与 与 文献 [2] 中 的 结果 
较为 一 致 。 
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图 10 火焰 中 心 紧 直 轴 线 的 温度 分 布 


Fig.10 Flame temperature profile along central vertical axis 
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